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Neutrino oscillations in the Earth matter may introduce peculiar modulations in the supernova 
(SN) neutrino spectra. The detection of this effect has been proposed as diagnostic tool for the 
neutrino mass hierarchy at "large" 1-3 leptonic mixing angle On- We perform an updated study on 
the observability of this effect at large next-generation underground detectors (i.e., 0.4 Mton water 
Cherenkov, 50 kton scintillation and 100 kton liquid Argon detectors) based on neutrino fluxes from 
state-of-the-art SN simulations. The average energies predicted by recent simulations are lower than 
previously expected and a tendency towards the equalization of the neutrino fluxes of the different 
flavors during the SN cooling phase appears. As a consequence, the detection of the Earth matter 
effect should be more challenging than expected from previous studies. In particular, for a typical 
galactic SN at 10 kpc, none of the proposed detectors may be able to detect the Earth modulation in 
the SN neutrino signal. Due to larger differences in neutrino (as opposed to antineutrino) energies 
and fluxes, in a 100 kton liquid Argon detector the Earth effect becomes observable at few kpc. 
Only for very close-by stars, that could evolve in SNe (like Betelgeuse at 0.2 kpc) all three detectors 
would clearly see the Earth signature. Finally, taking Icecube as co-detector together with a Mton 
water Cherenkov detector to monitor the Earth effect, it is not able to detect any sizable variation 
in the SN neutrino event rate for any galactic supernova. 

PACS numbers: Pq, 97.60.Bw 



I. INTRODUCTION 

Physical and astrophysical diagnostics via supernova (SN) neutrino detection in underground detectors represent 
a subject of intense investigation in astroparticle physics. In particular, a lot of attention has been devoted to 
possible signatures of the Mikheyev-Smirnov-Wolfenstein (MSW) matter effect [l| on the neutrino flavor evolution 
in supernovae (SNe) 0. Moreover in recent years, novel phenomena have been found to be important in the region 
close to the neutrinosphere, where the neutrino density is such that the neutrino-neutrmo interactions dominate the 
flavor evolution, causing large coherent conversions between different species (see, e.g., Q for a recent review). Flavor 
conversions in SNe could imprint peculiar signatures on the observable SN neutrino signal, sensitive to the neutrino 
mass and mixing, in particular to the unknown mass hierarchy [4] . 

Due to the uncertainties in the calculation of the primary SN neutrino spectra, it seems difficult to establish oscil- 
lation effects solely on the basis of differences in the observed spectra relative to theoretical expectations. Therefore, 
in the recent literature, the importance of relatively model-independent signatures has been emphasized, e.g. in asso- 
ciation with the prompt neutronization burst [B], with the rise time of the early neutrino signal or with matter 
effects associated to the shock- wave propagation at late times [7l-[Tl|. One unequivocal signature widely discussed 
would be the observation of the Earth matter effects 0, [H, , because they induce a characteristic energy-dependent 
modulation on the measured flux when SN neutrinos cross the Earth matter before being detected. Earth matter 
effects could be measured in a single detector, if it has enough energy resolution and statistics to track the wiggles in 
the observed energy spectrum induced by the neutrino oscillations in the Earth. In this context, a Fourier analysis of 
the SN neutrino signal has been proposed as a powerful tool to diagnose this modulation [13, Ell- Moreover, even if 
a single detector could not resolve these modulations, the comparison of the signals from shadowed and unshadowed 
detectors may allow one to diagnose the Earth effects [1^1 • 

The observability of the Earth matter effect largely depends on the neutrino average energies and on the flavor- 
dependent differences among the primar y sp ectra. At this regard, recent supernova simulations indicate lower average 
energies than previously expected (6l,[l7l[l8|. and a tendency towards the equalization of the neutrino fluxes of different 
flavors during the cooling phase, i.e. for post-bounce times t > 1 s |18,, JJj. Remarkably, this trend eases the agreement 
of theoretical expectations with SN1987A data [13: [HI- Motivated by this new input, we find worthwhile to reevaluate 
the detectability of the Earth matter effects in large detectors. In particular, we refer to three classes of detectors 
proposed for low-energy neutrino physics and astrophysics, viz. water Cherenkov (WC) detectors with fiducial masses 
of 0(400) kton @, [II], fiquid scintillation (SC) detectors with masses of 0(50) kton and fiquid Argon Time 
Projection Chambers (LAr TPC) with fiducial masses of O(IOO) kton [2^. These three detection techniques are the 
backbones of the European project LAGUNA (Large Apparati for Grand Unification and Neutrino Astrophysics) [1^ 
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and the LBNE (Long Baseline Neutrino Experiment) towards DUSEL (Deep Underground Science and Technology 
Laboratory) in the US [l^l- Moreover, the project of the Mton-class WC detector Hyper-Kamiokande is currently 
discussed in Japan [13] ■ In particular, WC and SC neutrino experiments are mostly sensitive to SN Vf. through the 
inverse beta decay process z?e + P n + e+. On the other hand, LAr TPC would have a high sensitivity to SN 
I'e's, through the charged current interactions of i/g with the Ar nuclei in the detector. Moreover, we also consider 
the Icecube ice Cherenkov detector [2I], taken as "co-detector" to monitor the Earth effect in comparison with the 
energy-integrated signal measured in a Mton WC detector [16]. 

For all these detectors, we find that no signature of the Earth matter effect is observable for a typical galactic SN 
at d ~ 10 kpc. Also in the more optimistic case of a close-by SN at d = 1 kpc, the chances to detect the Earth matter 
signature appear statistically weak in the antineutrino signal. Conversely a signal would show up in the v^, signal 
detectable at a LAr TPC. Finally, the Earth matter signal would be detectable with high significance in both neutrino 
and antineutrino signals for relatively close stars which might evolve into core-collapse SNe at unpredictable future 
times like Betelgeuse, Mira Ceti, and Antares (at d ^ 0.2 kpc), provided that the electronics of the detector will be 
able to cope with such a huge rate of events. These new results based on the state-of-art SN simulations seem to 
dramatically change the previous conclusions based on more "optimistic" choices for the primary SN neutrino fluxes. 

The plan of our work is as follow. In Sec. II we present the neutrino signal from SN hydrodynamical simulations 
we adopt as benchmark for our study. In Sec. Ill we characterize the SN neutrino flavor conversions and the Earth 
matter effect. In Sec. IV the features of our reference detectors (fiducial mass, cross sections and energy resolution) 
are described. In Sec. V we present our results on the detectability of the Earth effect in each of the reference 
detectors, and in the Icecube signal compared with a Mton WC detector. Comments on our results and conclusions 
are illustrated in Sec. VI. 



II. NUMERICAL MODELS FOR SUPERNOVA NEUTRINO EMISSION 



The un-oscillated distributions (with v = {v^, Ve, t^xi — J^p,, ^'r)}) are parametrized in energy and time as follows 



(1) 



where 4>u(t) is the energy-integrated neutrino number flux for each post-bounce time t, L^(t) the luminosity and 
{E^{t)) the mean neutrino energy. The function f^(E,t) is energy spectrum normalized such that J dE f^{E,t) = 1 
and parametrized as in [29j 
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The variables L^{t), {E^{t)) and a^(t) are extracted from SN simulations. 

In the following we adopt as input the recent SN simulations performed by the Garching group (G) -see [1] for 
references- and by the Basel/Darmstadt group (B/D) [l^. In Fig.[T]we show the initial luminosities Ly (upper panels), 
average energies {E^ (middle panels) and time-integrated energy-spectra in different time windows (lower panels) for 
the different flavors Ve, Ve and and for two different simulations: i) accretion phase of a 15 Mq progenitor from 
the Garching group (left panel); ii) accretion and cooling phases of a 18 Mq progenitor from the Basel/Darmstadt 
group (middle panels for the accretion and right panels for the cooling). 

In the former case, the accretion phase is clearly visible in the neutrino light-curve, lasting till t ~ 0.2 s and 
appearing as a pronounced hump in the electron (anti)neutrino luminosities. The relative time-integrated spectra for 
the accretion phase (with t G [0, 0.25] s) shows that flavor-dependent flux differences are large, with a robust hierarchy 
for the time-integrated neutrino number fluxes ($° = J dt (j),^{t)), namely > 3> Due to this situation, 
the accretion phase would represent the best time-window to detect effects of flavor conversions. 

The Basel/Darmstadt results show an accretion phase lasting till t ~ 0.4 s (middle panel). The main features of the 
neutrino signal are similar to what discussed in the Garching case. However, these simulations provide the late-time 
SN neutrino signal from the cooling phase too. This is represented in the right panels of Fig. [TJ During the cooling 
phase all neutrino flavors originate close to the neutron star surface, where the material is neutron rich, suppressing 
charged-current reactions for Fg. Therefore, the luminosities and the average energies of and become quite 
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FIG. 1: Time evolution of neutrino luminosities (upper panels), average energies (E^) (middle panels) and time-integrated 
energy spectra (lower panels) for (continuous curve), Ve (dashed curve), and Vx (dash-dotted curve). The left panels refer to 
the accretion phase of a Garching simulation for a 15.0 Mq progenitor, while the other ones to a Basel/Darmstadt simulation 
for a 18.0 AIq for the accretion (central panels) and cooling (right panels) phase. 



similar. As a result, the time-integrated neutrino fluxes in t e [1.0, 10.0] s are so close, that one expects to be difficult 
to measure oscillation effects in this case. 

In Table U we report the parameters of the time-integrated SN neutrino spectra for the Garching and the 
Basel/Darmstadt simulations distinguishing between the accretion (t < 0.25 s for G and t < 0.4 s for B/D) and 
the cooling phase. These parameters will be taken as benchmark for the evaluation of the Earth matter effect dis- 
cussed in the following. From this Table, it is evident that the neutrino average energies provided by the recent SN 
simulations are significantly lower than what assumed in the previous studies (see, e.g. Table 3 of [14]). 



III. FLAVOR CONVERSIONS AND EARTH MATTER EFFECT 



A. Neutrino mixing parameters 

We assume the three neutrino mass eigenstates separated by the following neutrino mass squared differences [3l| 

Aml^ - ml - 2 =. 2.35 x 10"^ eV^ , (4) 
Am% = ml - m\ = 7.54 x 10"^ eV^ . (5) 

The sign of Am^^jj^ is not determined yet. Therefore we consider both normal (NH, Amatm > 0) and inverted hierarchy 
(IH, Am^jjj^ < 0) scenarios. The mass eigenstates are related to the flavor eigenstates (i^e, ;/,-) by means of the 



4 



TABLE I: Spectral-fit parameters for the neutrino and antineutrino fluxes integrated over tlie accretion pfiase of the 15 Mq 
Garching progenitor (G) and both the accretion and cooling phases of the 18 Mq Basel/Darmstadt (B/D) progenitor. 
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mixing angles [31| 

sin^ 9i3 = 0.02 and sin^ 6^2 = 0.31 ; (6) 
the 2-3 mixing angle 023 is not relevant for our purposes since we are assuming equal i'^ and Vr fluxes. 



B. No Earth crossing 



The emitted SN neutrino flux is processed by self-induced and MSW oscillation eiTects during its propagation. The 
self-induced effects would take place within r 0(10"^) km from the neutrinosphere whereas the MSW transitions 
take place at larger radii, in the region r lO^-lO^ km. As the self-induced and MSW effects are widely separated 
in space, they can be considered independently of each other. 

We start considering the accretion phase. In NH, and for the spectral ordering of the accretion phase, no self-induced 
flavor conversion occu rs .321- 341. Instead, in IH potentially large self-induced effects could be expected for neutrinos 
and antineutrinos |32h34| . However, it has been shown using results both from Basel/Darmstadt group simulations 
and Garching group ones [ssUsol l that the trajectory-dependent multi-angle effects associated with the dense ordinary 
matter suppress collective oscillations in iron-core SNe, as expected when large trajectory-dependent phase dispersion 
are induced by the matter would suppress the collective phenomena [iflj . 

The neutrino fluxes can only undergo the traditional MSW conversions in SNe while passing through the outer 
layers of the star. Therefore, it is straightforward to com put e the Ve flux at the Earth in the different cases For 
"large" 6'i3, as recently measured by reactor experiments |4ll.l42j. in NH one finds 

Fp^ = cos2 6*12^^^"^ + sin^ 9i2F°^ and F^^ = . (7) 

Instead, in IH one gets 

F^^ = Fl and F,^ = sin^ e^F^^ + cos^ O.^F^^ . (8) 

Concerning the cooling phase, the characterization of the fiavor conversions is less straightforward. Indeed, it is 
expected that self-induced flavor oscillations would be no longer inhibited by the ordinary matter density and multiple 
spectral splits should occur [i^li^]. However, since the spectral differences among different flavors are not large (see 
Fig. [Ij, we numerically checked that the effect of the self-induced oscillations would produce a flavor equilibration 
among the different neutrino species (see also [13,14^). This would reduce the possibility to observe any signature of 
Earth matter effects. Moreover, it is expected that at t ^ 2 s the non-adiabatic effects associated with the matter 
turbulences in the supernova envelope would produce a smearing of the MSW flavor conversions [Tol . l46l |47| further 
reducing the spectral differences. In the following, since we are interested in the time-integrated signal during the 
cooling phase, we will neglect all these complicated effects and we will assume that also during the cooling phase the 
oscillated fluxes are described by Eqs. (I7])-(I8]). Our approach is conservative, since the experimental detectability of 
the Earth matter signature would be even more challenging than in our simplified scenario. 
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C. Earth crossing 

If the supernova is shadowed by the Earth for a detector [4^ , neutrinos will travel a certain distance through the 
Earth, and therefore will undergo Earth matter oscillations during their propagation. Since the neutrinos arrive at the 
Earth as mass eigenstates, the net effect of these oscillations can be written in terms of the conversion probabilities 
Pie — P{vi — > Ve)- For large ^13, one can obtain the neutrino fluxes at Earth as follows 

F®^{l-P^e)Fl+P2eFl (9) 

for neutrinos, and 

Fl^{l- P,e)Fl + P.eFl (10) 

for antineutrinos. Here P2e = P{i^2 — > i^e) and P2e = P{v2 i^e) while propagating through the Earth. The analytical 
expressions for P2e and P2e can be calculated for the approximate two-density model of the Earth [l^ . When neutrinos 
traverse a distance L through the mantle of the Earth, these quantities assume a very simple form [2I. [T2j: 

o ™ , ™ s n / ATO^sin26'i2 \ , , 

P2e = sin^ 019 + sm26VX sin(2eVX ~ 26*12) sin^ ^ L , 11) 

12 V 12 12) ^ 4£;sin26''i^ / ' ^ ' 



P2e = sin^ 9i2 + sin 26"^^ sin(2055 - 2^12) sin^ ( , Z.:^ L ] , (12) 



ArriQ sin 2^12 
iE sin 29',' 



where 0^ and ^"2 are the effective values of 612 in the Earth matter for neutrinos and antineutrinos respectively [49| . 
The Earth crossing induces a peculiar oscillatory signature in the neutrino energy spectrum. From Eqs. (|7|)- (jl0p . it 
results that the Earth matter effect would be present for antineutrinos in NH and for neutrinos in IH. 



D. Power spectrum of the Earth matter signal 



The typical event rate associated with an inverse beta decay process is oc E'^Ff^[E) since the cross section a oc E^. 
Whereas the distance between the peaks in the neutrino energy spectra associated with the Earth modulation, increases 
with energy, the peaks in the inverse energy spectrum are nearly equally spaced. In order to show this behavior, it is 
useful to recast the Ve and i>e fluxes in the presence of Earth matter as jlj| 
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where AF^ = F^^ is the energy-dependent difference of the primary neutrino spectra and = sin(26'5^ — 26*12); 

we also defined the mass squared difference in the Earth Atoq — Attiq sin 20i2 / sin 20^ and the "inverse energy 
variable" y = 12.5MeV/i?. Analogous definitions hold in the antineutrino sector in terms of the difference of the 
antineutrino spectra and of 0i2- The Earth matter effects on SN neutrino signal can thus be identified at a single 
detector through peaks in the Fourier transform of the inverse energy spectrum [T^ . The neutrino signal is observed 
as a discrete set of events. Following ^4], we define the power spectrum of the TV detected events as 
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This function is related to the continuous power spectrum G{k) of the events distribution function q{y) oc E'^F®{E), 
being 
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FIG. 2: Power spectrum G(fc) of the Earth matter signaL The left panel shows G{k) for a flux ratio varying in the range 
0.8 < $°^/$°^ < 1.5, while the right panel shows the impact of the change of the difference of the average energies in the range 
1 MeV < (E^J - < 3 MeV. 



Earth effects introduce peaks in the power spectrum at specific frequencies. In Fig. [21 for an illustrative purpose, we 
show the power spectrum G{k)oi the function q{y). In particular, we aim at discussing the dependence of G{k) on the 
ratio of the initial fluxes (left panels), and on the difi'erence among the neutrino average energies — (Ejy,) 

(right panels). We take as benchmark values the ones used in that refer to the 2003 Garching simulations (sof . 
In particular, we fix (E'pJ = 15 MeV and = = 3.0. In the left panel we assume [E^^) = 18 MeV and vary 
the ratio of the total fiuxes between 0.8 and 1.5. The two extremes cases are representative of the fiux 

ordering during the cooling and the accretion phase, respectively. The large peak in G{k) at low values of k is the 
dominant contribution due to the first two terms in Eq. ((T4| . The peak at fc ~ 90 corresponds to the oscillations in 
the Earth matter. Since Pg's have lower average energy, and due to the suppression at low-energy associated with the 
behavior of the cross section, the peak decreases significantly increasing the ratio <I>p__/<I>p^. In the right panel, 
instead, we fix — 0.8 and vary {EjjJ) between 15 and 18 MeV. As expected, since reducing the difference 

among the average energies E^AF^ becomes smaller, the peak in the power spectrum is consequently reduced. From 
this simple parametric study, we expect that for the SN models described in Sec. II the detection of the Earth peak 
would be more challenging than in the previous literature. In particular, since the average energies for different i/ 
species are lower and closer than previously assumed, this would produce a significant penalty factor in both the 
accretion and the cooling phase. We will present a quantitative estimation of this effect in the following Sections. 

Here instead we discuss the appearance of the double-peak feature shown in the left panel of Fig. [2] In Fig. [3] we 
refer to the uppermost and lowermost curves of the left panel of Fig.[5]as representative cases of a single-peaked (upper 
right panels) or a double-peaked (lower right panels) power spectrum, respectively. Instead, the left panels show the 
contributions of Pe and Dx to the observable signal E'^F®. It is worthwhile to notice that the power-spectrum G(fc) 
of the signal [Eq. ([16])] can be decomposed as 

G{k) = Gp, (fc) + Gp, (fc) + Gp,p, (fc) , (17) 

where the first contribution on the right-hand-side is associated with E^{1 — P2e)Fp , the second with E'^P2eF^ , and 
the third with the cross-correlation of the previous two. The three terms in the power-spectrum are peaked at very 
similar frequencies. Since Gp^p^(A;) can assume negative values, it can sometimes determine the appearance of the 
double-peak feature, especially when the two different terms become comparable as in the lower right panel. 
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FIG. 3: Left panels; Contributions of the different flavors to the observable Earth-modulated signal E'^F®^ corresponding to the 
uppermost and lowermost cases of left panel in Fig. [21 The continuous curve represents the total E^F® flux, the short-dashed 
one corresponds to _B^(f — P2e)F^^, while the long-dashed one to P2eF§^. Right panels: Different contributions to the 
power-spectrum G{k). The continuous curve represents the total power-spectrum G(k), the short-dashed one corresponds to 
Gi?^{k), the long-dashed one corresponds to Gn^ljt), while the dashed-dotted curve is for Gv^v^(k) (see text for details). 

IV. NEUTRINO DETECTION 

In this Section we describe the main aspects and ingredients of our calculations of supernova neutrino event rates. 
The oscillated SN neutrino fluxes at Earth F^^ must be convolved with the differential cross section erg for electron or 
positron production, as well as with the energy resolution function of the detector, and the efficiency e (that we 
assume equal to one above the energy threshold), in order to finally get observable event rates (Slj : 

Ne= F^®(7e®Re®£ ■ (18) 

We will now describe the main characteristics of four types of detectors we have used to calculate the signals in 
the presence of Earth matter effects, namely water Cherenkov detectors, scintillation detectors, liquid Argon Time 
Projection Chambers, and ice Cherenkov detector Icecube. 
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A. Water Cherenkov detectors 



In large WC detectors, the dominant channel for supernova neutrino detection is the inverse beta decay of electron 
antineutrinos^ 

De+p^n + e+ . (19) 

For this process, we take the differential cross section from [s^]- The total cross section grows approximatively as 
E^. We fold the differential cross sections for e+ production with a Gaussian energy resolution function of width A. 
The value of A is predominantly determined by the photocathode coverage of the detector. For our calculations we 
assume [Slj 

Awc/MeV = 0A7y/Ee/MeY , (20) 
where is the true positron energy. We assume as fiducial volume 400 kton [2^. 



B. Scintillation detectors 



In liquid SC detectors, the main channel for SN neutrino detection is the inverse beta decay of Pe's, the same as 
that in WC. However, here the positrons are detected through photons produced in the scintillation material. Since 
a larger number of photons can be produced in a SC detector, these have typically a much better energy resolution 
than the WC detectors. The energy resolution of the SC detectors is determined by the number of photo-electrons 
produced per MeV, which for this type of detectors is expected to be given by as good as [53| 

Asc/MeV = 0.07\/£^e/MeV . (21) 

Indeed, the energy resolution of a SC detector may be better by more than a factor of 6 than a WC. Since the 
Earth matter oscillations described in the previous section may get smeared out by the finite energy resolution of the 
detector, it is clear that the energy resolution plays a crucial role in the efficiency of detecting Earth effects. For our 
studies, we assume a fiducial mass of 50 kton [2J]. 



C. Liquid Argon Time Projection Chambers 



LAr TPC detectors would be particularly sensitive to SN electron neutrinos through their charged current interac- 
tions with Ar nuclei 

u^+^"Ar ^'^^K* +e- , (22) 

which proceed via the creation of an excited state of *'^K and its subsequent gamma decay. The Q-value for this 
inverse beta decay process is 1.505 MeV. The cross-section for this charged current reaction is taken from [54]. The 
one for leptons in LAr TPC has been calculated by the ICARUS collaboration which reports [55l] 

ALAr/MeV = 0.11 v/£^e/MeV + 0.02 £;e/MeV . (23) 
The fiducial volume for SN neutrino detection is taken to be 100 kton [56l |. 



D. Icecube 



A galactic SN v burst would be detectable in Icecube by a sudden, correlated increase in the photomultiplier count 
rate on a timescale on the order of 10 s (see Ref. [53| for a recent description). In its complete configuration and 
with its data acquisition system, IceCube has 5160 optical modules jsvl] and about 3 Mton effective detection volume. 



We will neglect the subleading neutrino interaction cliannels in tlie detectors, assuming tliat they can be separated at least on a statistical 
basis. 
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representing the largest current detectors for SN neutrinos. The SN neutrinos streaming through the antarctic ice 
interact mostly through + p — >■ n + e"*" reactions. While fine-grained detectors, like WC detectors, reconstruct 
individual neutrinos on an event-by-event basis, IceCube only picks up the average Cherenkov glow of the ice. The 
detection rate is given by [1, [16} 

/"OO 

7^p^ = / dE F^^ E,,i{E) aiE) , (24) 
Jo 

with E-cc\{E) being the energy released by a neutrino of energy E and cr{E) the inverse beta-decay cross section. 
All other detector parameters (angular acceptance range, average quantum efficiency, number of useful Cherenkov 
photons per deposited neutrino energy unit, average lifetime of Cherenkov photons, effective photo cathode detection 
area) have been fixed to the fiducial values adopted in [I6'], to which we address to for further details. 

V. DETECTING EARTH MATTER EFFECT 

A. Single detector 

We present our results about the detectability of the Earth effect. For our numerical calculations, we assume the 
mass-mixing parameters as in Eqs. (HHS]) and we will also assume that the path-length crossed by neutrinos in the 
Earth is L — 6000 km. Figure |4] shows the power-spectrum GN{k) [defined as in Eq. (|15p ] of the SN v signal in 
the WC (upper panels), SC (central panels) and LAr TPC detectors (lower panels). We discuss the results for three 
different SN distances from Earth (d = 10, 1,0.2 kpc). Since the Earth matter probability is time-independent, we 
consider as neutrino signal the time-integrated rate during the accretion phase taken from the 15 Mq Garching SN 
simulation (see Sec. 2 and Table U for the time-integrated flux parameters). We produce many realizations of the 
SN neutrino spectra at the Earth via Montecarlo simulations. The thin gray lines in Fig. 2] correspond to different 
realizations of G'Ar(fc) [Eq. P3|) ]. The thick black line corresponds to the power-spectrum G'Ar(fc) averaged over the 
different realizations. The light band around the average corresponds to the ±1ct level. 

Starting with a typical SN at d = 10 kpc (left panels), GAr(fc) ~ 1 for k ^ 40, for all the three detectors as expected 
in the absence of Earth modulation. Therefore, no peak in Gjv(fc) associated with the Earth effect is visible. For a 
SN at c? = 1 kpc, a peak around fc ~ 70 seems to emerge in the average GAr(fc) for the WC and for the LS. However, 
the power spectrum is compatible with 1 at la, i.e. with the expectations without the Earth crossing. Conversely, in 
the case of LAr TPC the peak in GAr(fc) is clearly visible at /c ~ 80. Indeed, this detector has the benefit of testing 
the Earth effect in the neutrino channel where the difference in the average energies/fluxes of v^^ and is larger than 
in the antineutrino sector. Therefore, the Earth signature is enhanced. We checked that the Earth effect starts to be 
visible in a LAr TPC at a few kpc for our benchmark SN neutrino fluxes. Finally in the right panels we show the 
power spectrum for the lucky case of a very close-by SN at d = 0.2 kpc. In this case, the peak in the power spectrum 
is clearly visible in the three different detectors. In particular, in the case of a liquid SC the superior energy resolution 
allows to see the double-peak structure in the power-spectrum discussed in Sec. 3. 

We repeated the same analysis with other eight SN models from the Garching simulations with progenitor mass 
between 12 and 40 Mq (see Fig. 1 in [6]) finding similar results to what shown before. Also for the accretion phase 
of the Basel/Darmstadt progenitors with 10.8 Mq and 18 Mq [T^, our results are similar. Finally, we calculated the 
power-spectrum in presence of the Earth matter also during the cooling phase for the Basel/Darmstadt simulations, 
characterizing v oscillations as described in Sec. 3. We find that, since during this phase the spectral differences 
among different flavors are very small, no peak in the power spectrum appears even in the optimistic case of a very 
close- by SN (results not shown here). 

B. Two detectors 

The Earth effect could produce a modification in the SN light-curve measured by Icecube. Therefore, together 
with a high-statistics Mton WC detector it could detect the Earth effect (if only one of the two detectors is shadowed) 
by the relative difference in the temporal signals [l6j . 

Figure [S] shows the ratio of the counting rate for IceCube and for a 400 kton WC detector, taking the input for 
the un-oscillated ly fluxes from the accretion phase of the 15 Mq Garching model. The counting rate in IceCube 
also includes the noise from the photomultipliers (280 Hz per optical module). We assume that the WC detector is 
un-shadowed. In the case that IceCube is shadowed we indicate the ratio of events with squares, while when Icecube 
is un-shadowed the ratio is plotted with circles. The left panel represents the ratio for a SN at 10 kpc, the central 
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FIG. 4: Power spectrum GN{k) of the Earth matter effect for a Galactic SN at d = 10 kpc (left panels), d = 1 kpc (central 
panels) and d = 0.2 kpc (right panels). The upper panels refer to a 400 kton WC detector, the middle panels to a 50 kton SC, 
and the lower panels to a 100 kton LAr TPC. The light curves corresponds to different MonteCarlo realizations of the power 
spectrum G'jv(fc), the thick curve to the average over the different realizations, and the band to the ±la variance level. 



at 1 kpc, and the right at 0.2 kpc, respectively. In all the cases considered in this section, we employ the actual 
time-dependent values of the parameters L^, (E^) and as obtained from 15 Mq Garching model. Note also that 
the shape of the ratio changes with the distance of the SN. In fact the measured signal in Icecube is the sum of 
a time-independent background rate (independent of the distance) and a true SN lightcurve whose normalization 
depends on the distance. 

The left panel for a SN at d = 10 kpc clearly shows that the ratio with and without Earth matter effects are nearly 
inseparable, since the statistical errors are larger than the difference. For a SN at d = 1 kpc the ratio in the shadowed 
and un-shadowed case still have too large statistical errors to be clearly separated. Conversely a.t d — 0.2 kpc the 
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FIG. 5: Events ratio in IceCube with respect to a 400 kton WC detector. The Mton detector is always considered un-shadowed 
by the Earth. The case when IceCube is shadowed is shown with the orange squares, whereas the case of Icecube un-shadowed 
plotted by black circles. The left panel represents the events ratio for a SN at 10 kpc, the central at 1 kpc, and the right at 0.2 
kpc, respectively. 

shadowed and un-shadowed cases are statistically separable. However, since the differences between the two curves 
are relatively small (maximum difference 1.7% at 40 ms) and the the initial neutrino fluxes are not known at this 
level, it is unlikely that the ratio of events could be used to diagnose the Earth effect, unless the comparison can be 
performed between two very similar detectors, with systematics canceling out. Otherwise, moderate variations in the 
distance of the SN events or the time-dependence used for the parameters {L^,, {E^) and oii,) may alone alter the 
"reference" curve with respect to which compared an eventual signal. 

VI. CONCLUSIONS 

Earth matter effect in supernova neutrinos would be in principle an interesting tool to diagnose the neutrino mass 
hierarchy at "large" 6*13. Motivated by the recent measurement of this angle |4l|, H^], and by the vivid discussion 
for future large underground neutrino detectors, we found worthwhile to reevaluate the chance to detect this effect 
in future neutrino experiments. In order to achieve a realistic characterization of this signature we used state-of- 
the-art SN simulations @, [l^ to describe the un-oscillated SN neutrino signal. The detection of the modulation 
in the neutrino spectra induced by the Earth crossing largely depends on the neutrino average energies and on the 
flavor-dependent differences between the primary spectra. At this regard, recent supernova simulations indicate lower 
average energies than previously expected [q. |l7l, |l8| , and a tendency towards the equalization of the neutrino fluxes 
of different flavors during the cooling phase [1^ . This makes the detection of the Earth matter effect more challenging 
than what assumed in previous works based on outdated SN simulations. 

In order to diagnose the modulation in the SN v energy spectrum induced by the Earth crossing, we perform a 
Fourier analysis of the neutrino signal in a 400 kton WC detector, in a 50 kton liquid SC and in a 100 kton LAr TPC. 
For all these detectors, we found that coming from a typical galactic SN at d = 10 kpc, no signature of the Earth 
matter effect is observable in the measured SN neutrino burst. Moreover, also in the more optimistic case of a close-by 
SN at d = 1 kpc, the chances to detect the Earth matter signature appear statistically weak in the antineutrino signal. 
Conversely, a signal would show up in the signal detectable at a LAr TPC. Only for relatively close stars which 
might evolve into core-collapse supernovae at unpredictable future times like Betelgeuse, Mira Ceti, and Antares (at 
d ^ 0.2 kpc), the Earth matter signal would be detectable with high significance in both neutrino and antineutrino 
signals. Finally, Icecube taken as co-detector to monitor the Earth effect together with a Mton WC detector is not 
capable to detect any sizeable variation in the SN neutrino event rate associated with the Earth matter effect for any 
galactic supernova. 

These new results based on the state-of-art SN simulations seem to dramatically change the previous conclusions 
based on more optimistic choices for the primary SN neutrino fluxes. In particular Mton WC detectors and large 
liquid scintillators would be able to observe the Earth matter signature only for very-close by (and rare) SNe, provided 
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that the electronics of the detector will be able to cope with huge rates of events. A 100 kton LAr TPC which starts 
to monitor the Earth signature from SNe at a distance of few kpc from the Earth, would statistically have ^ 10% of 
chance to see this signature from the next galactic supernova explosion [i^ . 

As a consequence of our finding, the possibility to determine the neutrino mass hierarchy with the next galactic SN 
neutrino burst requires to be rediscussed. Of course, a caveat is that, while current SN simulations have improved in 
many ways with respect to one or two decades ago, the results obtained should still be considered as indicative, and 
an empirical test would certainly be welcome. Turning the argument around, we can say that — barring an exceptional 
situation of a close-by SN — a positive detection of the Earth matter effect in future data would come as a surprise 
and probably invalidate the current models. A negative result most likely would turn into constrains in the flavor flux 
difference vs. average energy differences parameter space, thus indirectly testing these models. 

Of course, this should not discourage experimentalists to devote efforts to achieve a detailed measurement of the 
neutrino flux and spectra from a Galactic supernova; it would be "per se" a bonanza for testing the astrophysical 
models and the detailed understanding of the core-collapse SN mechanism. Also, matter effects in the Earth would 
still play an important role in determining the neutrino mass hierarchy, e.g. combining the NOvA long-baseline 
accelerator experiment ^5§\ with atmospheric neutrino data from the India-based Neutrino Observatory (INO) [s^ 
(see H). 

If Earth matter effect in SNe seems now not so promising as thought before, there are still other intriguing signatures 
in the SN neutrino signal that could give important information on the mass hierarchy. In particular, the early z/g 
ncutronization peak that could be detected in a Mton WC detector or in a large LAr TPC would provide a clear 
signature to extract the neutrino mass hierarchy Q . Also the early signal rise during the accretion phase measurable 
in all the detectors discussed in this work could encode an imprint of the neutrino mass ordering [6|. In conclusion, 
supernova neutrinos still represent a unique astrophysical probe of neutrino physics and astrophysics under extreme 
conditions. 
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